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Hadronic modes in the quark-gluon plasma and their spectral properties are discussed 
on the basis of the lattice QCD data. 
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' 1. Introduction 

I , Fate of light hadronic excitations above the critic al t emperature Tc of the QCD 

phase transition was first investigated 20 years agoPJ^ 
, In Ref.Q] whose title is "Fluctuation effects in hot quark matter: precursors of 

chiral transition at finite temperature" , the hadrinic modes above Tc were discussed 
in connection with the critical fluctuation associated with chiral symmetry restora- 
tion: The abstract says "■ • ■ there arise soft modes having a large strength and a 
' narrow width above the critieal temperature, which are analogous to the fluctuation 

of the order parameter in a superconductor above the critical points An example of 
such soft mode (ct-tt) at T > Tc was demonstrated using the Nambu-Jona-Lasinio 
modeP^ and was later called the "para-pion" From the modern point of view, It 
could be identified as a pre-formed qq resonance in the pseudo-gap phase where the 
chiral order parameter has a large fluctuation above Tc (see e.g. Refs. Ell- 
in Ref.El whose title is "Conjecture concerning the modes of excitation of the 
quark-gluon plasma" , the hadronic modes above Tc were discussed with a differ- 
ent perspective which is unique in non-Abelian gauge theories, i.e. the dynamical 
confinement of the plasma above Tc'. The abstract says " • ■ • the plasma exhibits 
confining features similar to that of the low-temperature hadronic phase. The con- 
fining features are manifest in the long-range, i.e., long-wavelength, low- frequency, 
modes of the plasma." 

Recently, studies of the hadronic modes in the time-like domain become possible 
in lattic e QC D simulations owing to the technique of the maximum entropy method 
(MEM).L^J2l It was then found that some hadronic modes (in the ss sector) survive 
even at T > Tc m quenched simulations (see Sec. 2))!^ However, its relation to the 
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above ideas (chiral pseudo-gap or dynamical confinement) is not yet clear. 

For heavy hadrons such as the charmonia, it has been considered that they 
dissolve soon after the deconfinement bec ause th e screened Coulomb potential above 
Tc is not st rong enou gh to hold the cc pair PIEl However, the quenched lattice QCD 
simulation^2EEljiave recently shown that resonances such as and rjc survive 
even up to T/Tc ^ 1.6 (see Sec. 0J. Physical origin of this phenomena is also an 
open question at the moment. 



2. QCD spectral functions and MEM 

The Matsubara correlation in a mixed representation is defined as Z3(r, p) = 
J dPx I?(t, x)e^*P '', where P(r, x) is the imaginary time two-point function. The 
spectral representation tells 

/ + 00 p"'^'^ 
A(^,p)d^ (0<r</3), (1) 
oo ^ ~r 

which is convergent for r 7^ in QCD. p) contains all the information of 
hadronic modes (mass, width etc) in the relevant channel. 

Reconstructing a continuous function A{ljj^ p) from discrete and finite lattice 
data i?(T, p) is a typical ill-posed problem. One way to avoid this difficulty is to 
introduce an ansatz for the spectral function with a few parameters. "'^'^ The max- 
imum e ntropy method (MEM) provides an alternative and much more powerful 
approacti I t one can obtain a unique A from the lattice data D without making 
a priori parameterizations. Also, statistical significance of the resultant A can be 
evaluated. 

In MEM, the most probable A given lattice data D is obtained by maximizing 
the conditional probability ^[AIL'] oc e"'^~^, where L is the standard likelihood 
function and S is the Shannon- Jaynes information entropy: 

(A{uj) 



A{uj) - ■m{uj) - A(w)log 



duj. (2) 



The statistical significance (error) of the resultant A is estimated by the second 
variation, [5 / 5Ay'P[A\D]. The default model m{ijj) in Eq.lj^J may be chosen so 
that the MEM errors become minimum. The final result is given by a weighted 
average over the parameter a as A{ijj, p) = / Aa{(-0, p)P[a|-Dm]da, where Aa{u), p) 
is obtained by minimizing P[^|Z)] for a fixed a. 

First successful application of MEM to the lattice QCD data at T = was 
carried out for the ground and excited mesons in Ref. i Also, basic concepts and 
techniques of MEM applied to lattice QCD are summarized in Ref. |H1 

3. Applications of MEM (T = 0) 

Shown in the left panel of Fig.^is a spectral image of the vector meson at rest (p = 
0) extracted from the quenched QCD data^The first (second) peak corresponds to 
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Fig. 1. Left panel: The dimensionless spectral function p(a;) = A(c.ij, 0)/(3w^) of the vector meson 
extracted from a 20^ X 24 lattice with /3 = 6.0. The hopping parameter is take to be k = 0.1557. The 
figure is adapted from Ref. 7. Right panel: The dimensionless spectral function p{u>) = A{uj, 0)/ui^ 
in the nucleon channel extracted from a 32^ lattice with /3 = 6.0 (a ~ 0.093 fm). The hopping 
parameter is take to be k = 0.1550. The horizontal axis denotes a dimensionless frequency Lua. 
This figure is adapted from Ref. 18. 



the ground (excited) vector meson. On the other hand, the highest peak corresponds 
to a bound state of Wilson doublers as first pointed out in Ref. El Shown in the 
right panel of Fig. ^ is t he s pectral function in the nucleon channel extracted from 
the quenched QCD data.li^ The first (second) peak corresponds to the nucleon and 
the Roper resonance, while the higher two peaks correspond to the bound states of 
Wilson doublers. 



4. Applications of MEM (T ^ 0) 

In Refs. |21 ^2 quenched QCD simulations with /3 = 7.0 on 32'^ x Nt anisotropic 
lattice were performed using naive plaquette gauge action and the standard Wilson 
quark action. The rcnormalized anisotropy is ^ = a^jar — 4.0 with a,- = ao./4 ~ 
9.75 X lO"*^ fm and — 1.25 fm. Employing the anisotropic lattice is necessary to 
keep the number of temporal data points as large as possible; N-t is taken to be 96, 
54, 46, 44, 42, 40 and 32 which correspond to T/T^ = 0.78, 1.38, 1.62, 1.70, 1.78, 
1.87 and 2.33, respectively. 

ss mesons at T 7^ 
Shown in FigI21 is the spectral functions of the ss mesons in the scalar (S), 
pseudo-scalar (PS), vector (V) and axial-vector (AV) channels for TjT^ = 1.38.1*^ 
The quark mass is chosen to reproduce the experimental ^ meson mass at T = 
approximately. One finds that spectral functions in all channels have degenerate 
peaks around 2.5 GeV at T/Tc = 1.38 (2.4 times the 0-meson mass at T = 0). 
Also, these peaks have been shown to disappear at higher temperatures e.g. at 
T/Tc = I.87P 

cc mesons at T ^ 
Shown in Figgis the case for charmonia.EI Quark masses are chosen so that we 
have rn^fi^ ~ 3.10GeV and m^^l ~ 3.03GeV at T = 0, which should be compared 
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Fig. 2. The dimensionless spectral function measured on an anisotropic lattice for the ss mesons 
at T/Tc = 1.38. The figure is taken from Ref. 9. 
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Fig. 3. Spectral functions measured on an anisotropic lattice for (a) the J/'^ channel {p{u)) = 
A{ll),0)/{3u)'^)) and for (b) the rjc channel {p{u>) = A(Lj,0)/iLj-^). Figures adapted from Ref. 12. 



with the experhiiental values, 3.10 GeV and 2.98 GeV, respectively. 

Sharp peaks for J/'$ and rjc are found even up to T ~ 1.62Tc as shown in Fig|3| 
with the peak position almost the same as their T = values. Similar observations 
on the charmonium bound states above Tc were obtained in other independent 
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our case, the sharp peaks disappear suddenly between l.GTc and 
1.7Tc. The width of the first peak in Fig|31partly reflects the unphysical broadening 
due to the statistics of the lattice data and partly reflects possible physical broad- 
ening at flnite T. The second and third peaks in Fig|31 are likely to be related to 
Wilson doublers. 

reliability of the MEM images 

Lack of accurate lattice data can cause fake peaks and/or fake smearing as 
demonstrated by the mock data in Ref. |H| Therefore the spectral functions obtained 
from MEM should pass some reliability tests. The first test is the error analysis of 
the peaks from the second variation, S'^P[A\D]/dA{LLj)6A{u!'). The second test is the 
change of the spectral functions under the variation of the number of data point 
-^data employed for the MEM analysis. The third test is a study of the finite volume 
effect in the spatial direction. The results in FigOlhave already passed the first and 
second tests. An attempt toward the third test is recently reported in Ref. El 

effects of dynamical quarks 

If dynamical quarks are included, there arise two competing effects to the disso- 
ciation rate of the hadronic modes observed in the quenched simulations: enhance- 
ment due to collisions with dynamical quarks, and suppression due to small critical 
temperature in fuU QCD (T^q""i<=h ^ 270 MeV vs. T^'"" - 170 MeV). The number 
ratio of the plasma constituents between full QCD {Nf = 2) and quenched QCD 
{Nf = 0) reads 

»q+g(ri"") ^ 16 + 21 / rj"" 

^^|-j,quench-| ^6 ^quench 

Therefore, the net thermal dissociation rate of the hadronic resonances may be even 
smaller in full QCD than that in quenched QCD for given T/Tc- In any case, it is 
truly necessary to carry out full QCD simulations to make definite conclusions. A 
preliminary study in Nf = 2 QCD is recently reported in Ref. 201 




5. Concluding remarks 

In this report, we have discussed the old idea of hadronic modes above Tc (chiral 
pseudo-gafPand dynamical confinemenlpl) and recent quenched lattice QCD results 
showing the existence of resonance structures above Tc in ss and cc channels. To 
unravel the true nature of such hadronic modes above Tc, it is important to extr act 
more information from lattice QCD data, such as the spatial wave functonsj^l 
spectral functions with finite spatial momentum,*^- role of had ronic modes to bulk 
plasma properties and q-q and o^otentials above T.P Further theoretical 
studies on the possible mechanism a I I to support hadronic states above Tc are 
also called for. 
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